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ABSTRACT 


EXTRAPOLATION  OF  INDUCED  TRANSIENT  SURFACE  CURRENTS  ON  A 
SCATTERER  TO  GENERAL  BAND-LIMITED  EXCITATION  CONDITIONS 

A technique  for  the  extraction  of  an  approximate  Singularity  Expan- 
sion description  of  the  electromagnetic  response  of  a scatterer  is  pre- 
sented. The  description  is  extracted  from  transient  waveforms  which  are 
spatial  samples  of  the  transient  current  on  the  scatterer  excited  by  a 
known  excitation.  The  method  can  use  any  reliable  Prony-type  processor 
for  exponential  analysis  of  the  transient  waveforms  so  long  as  it  pro- 
vides physically-meaningful  poles  from  the  waveforms.  The  fact  that 
Singularity  Expansion  poles  are  common  to  all  waveforms  on  the  object 
provides  a computational  advantage  over  a waveform-by-waveform  Prony- 
type  analysis.  A "consensus  pole  set"  procedure  is  described  to  exploit 
this  redundancy.  Some  issues  relating  to  the  completeness  of  a set  of 
transient  data  to  the  class  of  excitations  to  which  the  Singularity  Ex- 
pansion is  applied  is  discussed.  The  one-dimensional  example  of  a thin- 
wire  scatterer  is  used  to  test  the  practicability  of  the  method.  The 
ultimate  utility  of  the  method  is  demonstrated  by  way  of  expansion  of 
the  extracted  data  to  new  excitations.  Recommendations  are  made  toward 
an  implementation  of  the  method  using  measured  data. 
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CHAPTER  I 
INTRODUCTION 

1.1  SEM  Description  of  a System 

The  insight  that  the  general  transient  response  of  a scatterer 
resembles  a sum  of  exponentially  damped  sinusoids,  led  to  the  development 
of  the  Singularity  Expansion  Method  (SEM)  by  Baum  in  1971  [1,2].  The 
Singularity  Expansion  Method  has  since  become  a viable  procedure  for 
solving  for  the  transient  scattering  from  electromagnetic  systems  in  the 
terms  of  the  complex  resonant  frequencies  of  the  structure.  A collection 
of  SEM  parameters  for  a given  scatterer  include  the  complex  resonances 
along  with  the  associated  modal  current  distributions  and  normalization 
coefficients.  It  is  possible  to  expand  this  set  of  SEM  parameters  into  a 
transient  response  for  the  given  object  under  excitation  conditions 
limited  only  by  the  extent  of  the  parameter  set. 

To  date,  several  open  questions  concerning  the  SEM  representation 
have  lingered.  In  particular, 

-Of  the  several  coupling  coefficient  forms  postulated,  which  are 
valid? 

-Must  we,  in  general,  include  a contribution  to  the  transient  response 
which  results  from  an  entire  function  in  the  Laplace  transform 
domain? 

-How  do  the  issues  relating  to  the  entire  function  contributions  and 
coupling  coefficient  forms  interact? 

-What  means  exist  to  determine  the  SEM  description  for  a scatterer 
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whose  geometry  is  too  complex  to  admit  to  analytical  solution? 

The  first  three  of  these  questions  are  presently  under  investigation. 
The  present  work  is  an  effort  to  provide  a means  of  addressing  the  fourth. 

1.2  Prony-Type  Data  Analysis  as  Motivated  by  SEM 

Until  recently  the  only  method  by  which  the  SEM  parameters  were 
obtained  for  a structure  was  by  using  analytical /numerical  methods  such 
as  integro-differential  equations  or  boundary  value  techniques.  This  is 
practical  for  simple  object  geometries,  but  for  a complex  shaped  object 
such  as  an  aircraft,  these  approaches  are  impractical.  The  computational 
difficulty  for  the  method  is  overwhelming  for  complex  shapes.  Therefore, 
techniques  to  efficiently  extract  the  SEM  parameters  from  experimentally 
derived  data  are  sought. 

The  extraction  of  SEM  poles  and  residues  of  a system  from  its  tran- 
sient response  was  first  suggested  by  Mittra  and  Van  Blaricum  [3]  and 
extensively  elaborated  upon  by  them  [4,5].  The  development  of  a least- 
squares  Prony  procedure  by  Van  Blaricum  and  Mittra  [6]  gave  a method  to 
obtain  these  poles  and  residues  from  the  transient  response  of  the  scat- 
terer.  Though  the  least-squares  Prony  procedure  is  computationally  elegant 
and  reliable,  Dudley  [7]  has  shown  it  Introduces  a bias  in  the  estimate  of 
the  real  part  of  SEM  poles  when  noise  is  present  in  the  transient  data. 

An  alternative  technique  to  extract  the  SEM  poles  and  residues  from 
transient  responses  in  the  presence  of  noise  using  a Prony-related  proce- 
dure has  been  given  by  Lager,  et  al.  [8].  This  technique  is  the  so- 
called  "sliding-window  Prony"  procedure.  Work  by  others  to  find  a pole 
extractor  tolerant  of  noise  continues,  but  for  the  present  work,  only 
the  application  of  existing  pole  extraction  techniques  has  been  consid- 
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ered.  For  the  reader  interested  in  the  current  developments  of  alterna- 
tive methods.  Van  Blaricum  has  compiled  an  excellent  bibliography  [9]. 

A review  paper  on  Prony-type  analysis  as  applied  to  electromagnetic 
scattering  problems  appeared  recently  in  the  IEEE  Transactions  [10]. 

1.3  The  Contribution  of  Present  Work 

This  document  reports  the  results  of  an  effort  to  develop  a proce- 
dure to  extract  the  SEM  description  from  a transient  data.  Chapter  2 
presents  the  method  used  to  derive  the  SEM  description  from  spatially- 
sampled  transient  data.  Signal  processing  and  data  completeness  consid- 
erations are  also  discussed  there.  Chapter  3 reports  the  SEM  parameters 
which  have  been  extracted  from  transient  data  for  several  numerical  exam- 
ples. The  relationship  between  recoverable  SEM  data  and  the  frequency 
spectrum  of  the  excitation  are  discussed  there.  Some  results  obtained 
from  extrapolating  these  extracted  SEM  descriptions  to  new  excitation  con- 
ditions are  given  in  Chapter  4.  Chapter  5 presents  conclusions  drawn 
from  the  present  study  and  suggests  directions  which  future  work  in  this 
area  might  take. 
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CHAPTER  II 


DESCRIPTION  OF  THE  TRANSIENT  EXTRACTION  METHOD 


2.1  Description 


In  this  chapter  we  develop  a method  whereby  an  approximate  SEM  de- 
scription for  a given  scatterer  may  be  derived  from  spatial  samples  of 
transient  currents  excited  on  the  scatterer. 

Figure  1 pictures  a scatterer  with  unit  vectors  p..  symbolically 
indicating  current  sample  locations  and  orientation.  A known  incident 
field  EinC(r,t)  excites  a transient  current  J(r,t)  on  the  object.  The 
vector  r indicates  the  spatial  Independent  variable  and  t the  time 
variable.  Sampled  current  responses  at  the  sites  p^  are  scalar  functions 
of  time  of  the  form 


TjU)  * Pj  • J(r j ,t)  . (1) 

The  SEM  description  for  the  currents  on  the  object  can  be  written  as 

J(r,t)  - I Si  n<2)(t)  f^)  J±(r)  e8^  (2) 

where  the  S^'s  are  termed  "normalization  constants,"  f(s^)  is  the  Laplace 

transform  of  the  forcing  function,  J^(r)  is  a "natural  mode"  vector  and 

contains  the  spatial  distribution  information  associated  with  a given 

pole,  the  s^'s  are  the  complex  natural  frequencies  or  "poles",  and 

(2) 

(t)  are  the  "Class  2 coupling  coefficients"  and  are  a measure  of  how 
much  of  a mode  is  contained  in  the  current  response  for  a given  incident 


*.  * , *•  V- <■  . v * ; 


The  Class  2 coupling  coefficient  form  [1,11]  is  chosen  and  no  entire 
function  contribution  Is  Included  in  the  representation  (2).  It  has  been 
conjectured  that  the  Class  2 coupling  coefficient  yields  a zero  entire 
function  [12].  (It  is  pointed  out  that  the  present  work  does  not  depend 
on  the  validity  of  this  conjecture.)  After  the  leading  edge  of  the  inci- 
dent field  has  passed  over  the  structure,  i.e.  t > Tc>  the  Class  2 
coupling  coefficient  becomes  equal  to  the  Class  1 coupling  coefficient, 
which  is  a constant  with  respect  to  time.  Hence  we  observe  that  for 
t > T , (2)  admits  to  an  exact  temporal  representation  in  terms  of  decay- 
ing exponentials.  1 

If  the  sampling  process  indicated  in  (1)  is  applied  to  the  SEM  repre- 
sentation for  the  current  given  by  (2)  the  result  is 

rj(t)  - I Bi  n[2\t)  f(s±)  Pj  • e8!*"  (3) 

For  time  after  the  leading  edge  of  the  incident  wave  has  passed  over 

the  structure,  i.e.  t > T , the  Class  2 coefficients  are  equal  to  the 

c 

Class  1 coefficients  and  the  sampled  currents  take  the  form  of  a sum  of 
decaying  exponentials  with  constant  coefficients: 

n^(t)  - const.  - t > Tc,  (4) 


and 


rj(t)  ■ l Au 


.•i* 


(5) 


where 


Aij  ■ h "l1'  *<•,> 


W- 


(6) 


From  here  forward  we  adopt  the  common  SEM/Prony  vernacular  and  term  the 


A^' 8 "residues"  and  the  s^'s  "poles".  This  terminology  is  based  on  the 
role  of  these  constants  in  the  Laplace  transform  domain.  Observe  r ^ (t) 
are  sampled  currents  which  are  real  functions.  Therefore,  A^'s  and  s^s 
must  occur  in  complex  conjugate  pairs  or  lie  on  the  real  axis. 

The  key  to  the  extraction  procedure  is  the  recognition  that  a func- 
tion of  the  form  (5)  can  be  decomposed  into  its  pole/residue  components 
by  means  of  Prony-type  analysis  [6].  Suppose  we  have  available  noise- 
corrupted  observations  r°(t)  of  the  r^ (t)  signals.  A Prony-type  analysis 
approximates  the  corrupt  signals  as 

r?(t)  * I A; 


O/^v  ..  Sijt 


ij 


(7) 


where  the  A.,  and  s are  the  output  of  a Prony-type  process  applied  to 


ij 


ij 


r^(t).  The  superscript  "o"  is  used  to  indicate  "observed"  values. 

In  order  to  determine  the  SEM  description  from  these  data,  we  must 
match  A°j  and  s°^  with  A^  and  si  from  equation  (5).  From  equation  (2) 
we  see  the  poles  of  a system  are  independent  of  sample  location  j . If 
the  observed  data  and  the  processing  thereof  were  perfect  this  result 
would  be  evidenced  by  s°j  ■ s°^  for  all  j,k.  Because  of  the  noise  in 
both  the  data  and  the  computer  processing  this  result  is  not  achievable. 
The  next  section  describes  a process  which  provides  a single  concensus 
extracted  pole  set  for  each  set  of  spatial  samples.  This  concensus  pole 
set  represents  the  first  of  the  extracted  SEM  quantities. 

The  remaining  SEM  quantities  are  imbedded  in  the  A^'s.  From  (6) 
we  write  in  terms  of  "observed"  SEM  quantities 


6°  n'1'  fCj)  Pj-jJftj)  - 


(8) 


The  observation  that  the  first  three  terms  on  the  left-hand  side  are 


8- 


constant  with  respect  to  the  spatial  variables  leads  to  the  conclusion 


that  the  left-hand  side  Is  a scaled  natural  mode.  Thus  we  construct 


tI,- 


Vr)  “ Ci  1 UilPl’  A12P2 AimPm} 


where  7 Is  an  Interpolation  operator  and  a scale  factor  which  Is  chosen 


to  yield  a sensible  normalization  of  J*(r) . 


With  this  approximate  knowledge  of  an  approximation  to  the 


coupling  coefficient  is 


n(1)o  - //  I (r,s)  • jJCr)  dS, 

object 


where  EQ(r,s)  is  the  propagation  factor  portion  of  the  excitation  field, 


i.e.  we  factor  the  total  excitation  as 


E(r,s)  - f (s^  E (r,s) , 


Finally,  the  constants  B°  may  be  determined  from  (8)  since  all  other 


factors  in  the  expression  are  known.  Viz. 


1 -41)0  ifSj)  fj  • ' 


It  should  be  pointed  out  that  the  index  j appears  in  the  right-hand  side 


of  (12).  However,  p^  • J°(r^)  is  proportional  to  so  the  ratio  in  (12) 


is  constant  with  respect  to  j . 


In  summary  the  key  steps  in  the  extraction  are  as  follows: 


1.  Apply  a Prony-type  algorithm  to  spatially-sampled  transient 


currents . 


2.  Construct  an  approximate  natural  mode  from  the  spatial  variation 


of  the  residues  obtained  for  a single  pole  per  (8)  and  (9). 
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3.  With  the  coupling  coefficient  computation  as  an  intermediate 
step,  compute  the  normalization  constants  8^. 

The  resulting  data  within  the  limitations  of  completeness  outlined  in  the 
following  section  may  be  reexpanded  to  provide  transient  currents  under 
new  excitation  conditions.  This  reexpansion  can  be  carried  out  formally 
with  Class  1 or  Class  2 coupling  coefficients  since  the  modes  and  normal- 
ization constants  are  known  explicitly. 

2.2  Signal  Processing  Considerations 

Most  data  analyzed  by  this  process  is  likely  to  be  plagued  by  noise 
whether  obtained  numerically  or  experimentally.  Dudley  [7]  has  shown 
when  the  least-square  Prony  procedure  is  applied  to  data  containing  noise, 
the  resulting  poles  are  biased  in  the  estimate  of  the  real  part.  There- 
fore, if  a Prony-type  technique  is  to  be  used  as  a practical  tool  it  must 
be  tolerant  of  noise  in  the  data  it  analyzes. 

We  know  that  a given  object  possesses  poles  which  are  intrinsic  to 
it  alone.  However,  the  pole  sets  computed  from  the  many  sampled  current 
waveform  observations  along  the  structure  are  not  identical.  Therefore, 
a technique  is  used  to  form  the  best  observed  pole  set  for  the  structure 
which  we  term  the  "consensus  pole  set." 

To  form  a consensus  pole  set  we  follow  the  steps  below: 

1.  Use  a pole  extractor  to  obtain  pole  sets  from  each  observed 
waveform; 

2.  Display  these  pole  sets  graphically,  and  group  only  the  poles 
which  tend  to  form  a cluster,  as  Illustrated  in  Figure  2; 

3.  Take  the  average  value  of  the  poles  grouped; 

4.  These  spatially  averaged  poles  and  only  these  averaged  poles 
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constitute  the  consensus  pole  set.  Any  other  poles  present  are 
presumed  to  be  "curve  fitting"  poles  [6]  and  are  discarded;  and 
5.  Use  this  consensus  pole  set  alone  to  obtain  corresponding  resi- 
dues for  each  sampled  waveform  and  the  remaining  SEM  description. 
This  procedure  is  essentially  the  "cluster  averaging"  method  used  by 
Lager,  et  al.  (8)  in  their  sliding-window  Prony  technique.  It  is  applied 
here  to  exploit  spatial  rather  than  temporal  redundancy  in  the  data. 

.The  pole  extractors  used  in  this  work  were  the  least-squares  Prony 
procedure  after  Van  Blaricum  and  Mittra  [6],  and  the  sliding-window  Prony 
procedure  after  Lager,  et  al.  [8].  Because  of  its  computational  effi- 
ciency, the  least-square  Prony  procedure  was  used  whenever  the  data 
exhibited  a low  level  of  noise.  A set  of  poles  was  obtained  for  each 
waveform.  These  pole  sets  were  associated  as  described  above  to  yield 
the  consensus  pole  set.  However,  when  noise  levels  were  encountered  with 
magnitudes  severe  enough  to  make  clusters  unidentifiable,  the  sliding- 
window  Prony  procedure  was  used.  This  occurs,  for  example,  in  numerical 
data  when  the  wire  is  locally-driven  near  its  end  where  coupling  is  weak. 

In  summary,  a set  of  poles  for  each  sampled  waveform  is  obtained  by 
an  appropriate  pole  extraction  technique.  The  poles  are  graphed,  those 
which  aggregate  in  clusters  are  grouped  and  averaged.  The  averaged  poles 
and  only  these  poles  constitute  the  consensus  pole  set.  The  residues  are 
calculated  for  an  expansion  in  terms  of  the  consensus  pole  set  alone. 

2.3  Data  Completeness  Considerations 

The  general  SEM  description  for  a system  consists  of  an  infinite 


collection  of  data.  However,  when  applied  to  measured  or  computationally 
derived  results  only  a finite  collection  of  data  can  be  obtained.  There- 
fore, the  question  of  data  completeness  arises. 
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In  general,  one  would  likely  want  to  acquire  sufficient  data  to 
acquire  spatlally-complete  data  within  a predefined  bandwidth. 

There  are  two  mechanisms  which  limit  the  extent  of  the  SEM  data 
derivable  from  a given  collection  of  transient  data:  the  extent  of  the 
frequency  spectrum  of  the  time  history  of  the  excitation;  and  the  coupling 
between  the  spatial  form  of  the  excitation  and  individual  modes. 

The  band-limiting  of  the  frequency  spectrum  limits  the  resonances 
which  are  excited.  Hence,  if  one  wishes  to  recover  high-order  poles,  he 
must  use  an  excitation  which  is  sufficiently  broadband  to  excite  the 
resonances  desired  to  an  amplitude  well  above  the  noise  level  of  the 
data. 

The  spatial  placement  and/or  distribution  of  the  excitation  can  be 
such  that  it  has  a small  or  zero  coupling  coefficient  for  certain  modes. 
For  example.  Figure  3 illustrates  a center-driven  cylinder  and  the  first 
antisymmetric  mode  on  a cylindrical  structure.  The  generator  is  located 
at  a null  of  the  mode.  Thus,  in  the  coupling  coefficient  integral  of 
(10),  this  excitation/mode  combination  yields  a zero  coupling  coefficient. 
Even  if  the  generator  were  slightly  off  center,  the  coupling  coefficient, 
though  non-zero,  would  be  small.  The  associated  resonance  would  be  only 
weakly  excited.  Hence  it  could  be  lost  in  noise. 

The  end  in  obtaining  the  SEM  description  is  to  be  able  to  extra- 
polate to  new  excitations.  However,  if  the  description  is  incomplete, 
either  in  the  spectral  or  the  spatial  sense,  then  extrapolation  cannot  be 
carried  out  for  a general  excitation.  In  the  spectral  sense,  the  extra- 
polation cannot  be  extended  beyond  the  band  limitation  of  the  description. 
In  the  spatial  sense,  the  SEM  description  cannot  be  expanded  to  situations 
which  properly  couple  data  not  present.  For  example,  if  an  even  excita- 
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Cion  were  used  Co  produce  Che  daCa  on  which  Che  exCraccion  process 
operaCes,  Chen  Che  SEM  descripCion  will  con Cain  only  even  modes  and  can 
be  excrapolaced  Co  even  exciCaCions. 


CHAPTER  III 


EXTRACTED  SEM  RESULTS 


3.1  Introduction 

This  Chapter  reports  the  SEM  results  extracted  from  several  different 
collections  of  transient  data  using  the  method  described  in  the  pre- 
ceding Chapter. 

The  transient  data  were  obtained  using  th>;  time  domain  computer  code 
TWTD  [13].  The  structure  modeled  by  TWTD  was  a thin  cylindrical  scat- 
terer.  This  structure  was  chosen  because  of  the  availability  of  compara- 
tive SEM  results  obtained  by  Tesche  [14]  and  because  of  the  simplicity 
incumbent  with,  the  one-dimensional  structure. 

A 1-meter  dipole  antenna  with  a half-length  to  radius  ratio  of  100 
was  modeled  by  TWTD  using  49  spatial  segments.  The  excitation  was  defined 
to  be  a local  generator  occupying  a single  zone.  Figure  4a  illustrates 
the  segmentation,  and  shows  the  numbering  scheme  used. 

For  excitation,  a Gaussian  waveform  in  time  was  chosen,  having  the 

form 


yt> 


/t-^m  \ 2 
- e (t72  } 


where  Tm  is  a time  shift  chosen  to  produce  quasi-casuality,  and  Tw,  termed 
the  waist  of  the  Gaussian,  is  used  to  adjust  the  spectrum  of  the  excita- 
tion in  the  frequency  domain.  As  evidenced  by  Figure  4b,  the  waist  of 
the  Gaussian  is  defined  as  the  width  of  the  Gaussian  at  36.8  percent  of 
its  peak  value,  or  the  1/e  point. 


Figure  4.  (a).  TWTD  model  of  the  wire  scatterer  with  the  generator  at 

the  center  segment:  zone  25.  The  order  of  segmentation  is  also  shown, 
(b).  Model  of  the  Gaussian  generator  used  for  excitation,  showing  the 
"waist"  parameter. 


Figure  5.  The  spectrum  of  the  Gaussian  forcing  function,  with  a waist  of 
0.6158  nanoseconds,  shown  with  the  poles  extracted  from  the  transient 
response  it  generated:  Example  1.  Note  the  scattering  of  the  poles 
above  pole  7. 
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are  identifiable  through  pole  17.  This  is  commensurate  with  the  spectrum 
of  the  excitation. 

3.3  Extracted  SEM  Parameters  - Even  Mode  Case 

The  extraction  technique  described  in  Chapter  II  has  been  applied 
for  the  cases  described  in  the  previous  section.  The  identifiable  poles 
were  grouped  and  averaged  to  form  consensus  pole  sets.  Residues  were 
extracted  according  to  the  least-square  Prony  procedure  using  the  consen- 
sus pole  set;  i.e.  the  "curve  fitting  poles"  were  discarded  and  cluster 
average  values  used.  The  data  which  follow  compared  with  those  given  by 

Tesche  [14].  Tesche's  data  comprise  only  modes  1 through  12.  Thus  no 

/ 

comparisons  are  possible  for  modes  13  through  17  in  the  second  case 
studied. 

The  consensus  poles  and  the  poles  reported  by  Tesche  are  compared  in 
Table  1.  The  nine  poles  determined  in  Example  1 agree  well  with  Tesche's. 
Notice  that  we  included  the  ninth  cluster  even  though  there  is  some  ques- 
tion as  to  whether  this  cluster  is  "identifiable."  Subsequent  results, 
however,  evidence  poor  quality  in  the  mode  corresponding  to  this  pole. 

The  poles  from  Example  2 compare  well  to  those  reported  by  Tesche, 
and  exhibit  a consistent  trend  for  the  poles  beyond.  However,  no  compar- 
ison is  available  for  these  higher  frequency  poles  so  the  exact  accuracy 
is  not  known. 

Figures  7 through  12  are  comparisons  of  the  extracted  modes  with 
those  of  Tesche.  Every  mode  obtained  for  Example  2 is  given,  while  only 
selected  modes  are  given  for  Example  1. 

Figure  7 displays  mode  1 and  7,  corresponding  to  poles  number  1 & 7 
respectively  for  Example  1.  All  of  the  odd-numbered  modes  of  Example  1 
are  essentially  indistinguishable  from  those  of  Tesche.  Figure  8 
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Consensus  Poles  and  Those  Reported  by  Tesche 


Pole 

Index 

Tesche 

sL/c 

Real  Imag 

POLE  VALUES 

Example  1 
sL/c 

Real  Imag 

Example  2 
sL/c 

Real  Imag 

1 

-0.082 

0.926 

-0.082  ■ 

-0.918 

-0.082 

0.918 

3 

-0.147 

2.874 

-0.148 

2.878 

-0.148 

2.878 

5 

-0.188 

4.835 

-0.189 

4.844 

-0.189 

4.844 

7 

-0.220 

6.800 

-0.222 

6.789 

-0.221 

6.796 

9 

-0.247 

8.767 

-0.256 

8.710 

-0.248 

8.717 

11 

-0.270 

10.733 

. 

-0.278 

10.600 

• 13 

-0.317 

12.392 

15 

-0.372 

14.101 

17 

-0.435 

15.714 
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compares  the  extracted  mode  9 with  Tesche's  corresponding  mode.  The  com- 
parison is  seen  to  be  unsatisfactory,  particularly  in  the  imaginary  part. 
The  excitation  spectrum  is  too  small  at  this  pole  to  provide  a recover- 
able mode.  Recall  that  the  spectrum  significantly  excited  only  through 
the  seventh  mode  and  it  is  questionable  as  whether  the  9th  pole  should 
have  been  included  in  the  consensus  pole  set  in  the  first  place. 

The  modes  of  Example  2 are  shown  in  Figures  9 through  12.  Modes  1 
through  11  agree  favorably  with  Tesche's,  but  since  he  reported  only 
these  eleven  modes,  comparison  of  higher  modes  is  not  possible.  However, 
the  higher  modes  show  the  features  one  would  expect. 

Table  2 gives  a comparison  of  Tesche's  normalization  constants  with 
those  of  both  the  preceding  Examples.  We  see  from  the  Table  that  the 
normalization  constants  from  Example  1 compare  well  with  Tesche's  through 
pole  7.  The  results  of  Example  2 agree  with  Tesche's  results  favorably 
for  all  the  poles  reported.  The  departure  between  the  two,  though  satis- 
factory,  cannot  be  overlooked  in  constants  9 and  11.  The  constants  seem 
to  grow  too  rapidly  for  poles  13  and  beyond.  The  precise  cause  for  this 
phenomenon  has  not  been  investigated. 

3.4  Extracted  SEM  Parameters  - General  Case 
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TABLE  2 

Transient  Extracted  Normalization  Constants  and  Those 
Reported  by  Tesche 


Pole 

Index 


Tesche 

L/c 

Mag  Ang 


NORMALIZATION  CONSTANTS 
Example  1 


Example  2 


15.59 

16.45 

15.86 

17.25 

18.83 

21.74 


3.68 

12.33 

3.68 

12.33 

4.75 

5.24 

4.75 

5.15 

5.53 

0.94 

5.62 

0.83 

6.18 

2.77 

6.19 

5.45 

5.03 

32.66 

7.27 

14.13 

7.87 

27.20 

9.84 

30.85 

11.54 

30.11 

13.18 

23.22 
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least-squares  Prony  procedure  was  first  applied,  but  consensus  clusters 
could  not  be  identified.)  Example  3 contains  poles  1,  2,  4,  5,  7,  9, 
and  10,  since  the  other  poles  are  not  coupled  spatially  for  this  excita- 
tion. In  Example  4,  poles  1,  2,  3,  5,  6,  7,  9,  and  10  are  coupled.  Com- 
paring pole  values  of  Examples  3 with  Tesche's  poles  we  find  good  agree- 
ment for  every  pole.  However,  some  small  differences  are  noted  in  the 
comparisons  of  Example  4. 

The  natural  modes  for  pole  numbers  1,  2,  and  5 are  given  in  Figures 
13  through  15,  for  each  of  these  Examples.  Extracted  modes  beyond  mode  4 
represent  the  poorest  quality  results  obtained  in  this  study.  The  essen- 
tial oscillatory  features  still  appear  in  the  modes,  but  the  zero  crossings 
are  shifted  to  a small  degree  and  errors  in  peak  amplitudes  on  the  order 
of  twenty  percent  are  present.  Precise  modal  symmetry  or  antisymmetry  is 
generally  lost  by  way  of  these  two  types  of  errors.  Observing  mode  1 of 
both  Examples  shown  in  Figure  13,  we  see  that  the  mode  of  Example.  3 is 
more  nearly  correct  than  that  of  Example  4.  We  conjecture  this  is  due  to 
the  location  of  the  generator  for  Example  3;  i.e.,  mode  1 of  Example  3 is 
more  strongly  coupled  than  mode  1 of  Example  4 and  hence  better  data  is 
derivable.  Figure  14  displays  mode  2 for  both  Examples.  For  this  mode, 
the  Example  4 result  is  more  strongly  coupled.  Evidence  of  this  is  seen 
in  the  more  stable  imaginary  part  and  proper  zero  crossings.  Figure  15 
shows  mode  5 for  both  Examples.  Here  we  see  results  for  both  Examples 
which  are  not  particularly  gratifying.  Example  4 is  seen  to  possess  more 
nearly  uniform  amplitude  peaks.  We  again  conjecture  this  is  due  to  the 
slightly  stronger  coupling  of  this  mode  in  Example  4.  Reexpanded  results 
in  the  next  Chapter  show,  however,  that  these  errors  do  not  strongly  effect 
extrapolated  results  so  long  as  the  spectrum  for  the  reexpanslon  excita- 
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TABLE  3 


Consensus  Poles  and  Those  Reported  by  Tesche 


Pole 

Index 

Tesche 

sL/c 

Real  Imag 

POLE  VALUE 

Example  3 
sL/c 

Real  Imag 

Example  4 
sL/c 

Real  Imag 

1 

-0.082 

0.926 

-0.078 

0.934 

-0.092 

0.922 

2 

-0.120 

1.897 

-0.124 

1.862 

-0.159 

1.916 

3 

-0.147 

2.874 

* 

-0.139 

2.789 

4 

-0.169 

3.854 

-0.168 

3.918 

* 

5 

-0.188 

4.835 

-0.189 

4.845 

-0.168 

4.807 

6 

-0.205 

5.817 

* 

-0.259 

5.853 

7 

-0.220 

6.800 

-0.230 

6.812 

-0.221 

6.797 

9 

-0.247 

8.767 

-0.248 

8.718 

-0.265 

8.926 

10 

-0.260 

9.752 

-0.272 

9.598 

-0.284 

9.668 
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tion  does  not  differ  drastically  from  that  of  the  source  data. 

Table  4 contains  the  magnitudes  and  phases  of  Tesche's  normalization 
constants  and  those  of  both  Examples.  Comparing  the  constants  of  these 
examples  to  Tesche’s  we  find  relatively  good  agreement  for  constants 
associated  with  pole  number  1.  Comparing  the  constants  of  pole  number  2, 
we  see  the  value  for  Example  3 is  much  closer  to  Tesche's  value  than 
Example  4.  These  errors  may  be  traced  to  the  larger  pole  errors  present 
in  the  Example  4 result.  For  example,  the  real  part  pole  value  for  the 
second  pole  is  in  error  by  approximately  32.5  percent,  resulting  in  a 
normalization  constant  whose  magnitude  is  in  error  by  38.9  percent. 

The  errors  in  the  normalization  constants  beyond  the  fifth  one  are 


quite  large  - several  hundred  percent  at  times.  This  is  due  largely  to 
the  fact  that  they  are  derived  from  a localized  excitation  and  that  the 
local  errors  in  the  modes  are  often  large  in  these  examples.  The  use  of 
several  generator  locations  in  a given  extraction  procedure  and  disturbed 
excitation  is  expected  to  control  this  source  of  error. 

The  results  associated  with  pole  5 in  Example  3 demonstrate  another 
mechanism  for  the  intrusion  of  error.  The  extracted  pole  value  Itself  is 
satisfactory.  However,  the  generator  location  is  such  that  the  mode  con- 
tributed only  weakly  in  the  total  response.  Observation  of  Figure  15 
shows  how  the  mode  for  this  pole  is  corrupted  because  of  the  nearness  of 
its  contribution  to  the  noise  level.  This  error  in  the  mode,  in  turn, 
propagates  into  the  evaluation  of  the  normalization  constant. 
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TABLE  4 

Transient  Extracted  Normalization  Constants  and  Those 
Reported  by  Tesche 


NORMALIZATION  CONSTANTS 

Pole  Tesche  Example  3 Example  4 

Index  L/c  L/c  L/c 

Mag  Angle  Mag  Angle  Mag  Angle 

1 3.72  ] 

2 4.21  ] 

5 5.32  ] 
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CHAPTER  IV 

EXTRAPOLATED  RESULTS 

4.1  Introduction 

The  utility  of  extracted  SEM  results  lies  In  the  ability  to  apply  it 
to  circumstances  other  than  those  through  which  the  data  were  derived. 

This  Chapter  presents  results  obtained  when  the  data  described  in  the 
previous  Chapter  were  used  to  characterize  the  thin-wire  scatterer's 
response  to  several  plane  wave  excitations. 

For  this  extrapolation  study,  we  created  a complete  SEM  data  base 
from  transient  extracted  parameters.  The  data  base  was  constructed  from 
transient  data  excited  by  a Gaussian  generator  with  a waist  of  0.3077 
nanosecond.  The  quantities  associated  with  poles  1,  3,  5,  7,  and  9 were 
derived  from  data  excited  by  a generator  located  at  the  center  of  the  wire, 
those  associated  with  poles  2,  4,  and  10  from  data  excited  by  a generator 
at  0.367  L,  and  those  associated  with  pole  6 from  the  data  excited  by  a 
generator  at  0.245  L.  Note  that  this  data  base  does  not  include  pole  8. 

The  cases  presented  here  do  not  require  it  in  their  representation.  It 
may  be  derived  readily  from  appropriate  transient  data.  It  was  omitted 
for  convenience  here. 

4.2  Gaussian  Plane  Wave  Excitation 

The  transient  extracted  SEM  data  base  was  reexpanded  using  a plane 
wave  with  a Gaussian  time  history  as  the  form  of  excitation,  for  0,  30, 
and  60  degrees  incidence,  and  the  current  was  observed  for  each  of  these 
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angles  at  observation  segments  0.245  the  distance,  0.367  the  distance,  and 
half  the  distance  from  the  end  of  the  antenna.  Two  comparisons  were  used 
to  check  the  validity  of  the  transient  extracted  SEM  reexpansion.  The 
first  is  computed  data  for  a direct  time  domain  solution  to  the  integral 
equation  for  the  wire  using  the  TWTD  program.  The  second  is  a reexpansion 
in  terms  of  SEM  parameters  obtained  by  a method  of  moments  computation 
similar  to  that  of  Tesche. 

Figure  16  displays  transient  currents  computed  by  each  of  these  means 
at  0 degrees  incidence  for  the  three  observation  points.  Here  we  observe 
almost  perfect  agreement  among  the  three,  except  in  the  early  time.  The 
reason  for  the  disagreement  in  the  early  time  is  that  only  the  pole  con- 
tributions due  to  poles  of  the  scatterer  are  included  in  the  SEM  expansions. 
The  Gaussian  excitation  requires  the  inclusion  of  an  entire  function  con- 
tribution in  the  inversion  process.  These  data  were  generated  primarily 
to  observe  the  late  time  behavior  of  the  extrapolated  results.  Conse- 
quently, this  contribution  was  not  included.  A subsequent  example  for  a 
different  generator  waveform  allows  early-time  comparisons  as  well. 

Figure  17  compares  the  currents  from  the  three  computations  for  30 
degrees  incidence  at  the  three  observation  points.  In  this  Figure  we 
again  see  good  agreement  in  the  late  time.  This  comparison  also  manifests 
departures  of  the  SEM  expansion  from  the  TWTD  computation  in  the  early 
time. 

The  currents  excited  by  a plane  wave  incidence  of  60  degrees  are  shown 
in  Figure  18.  The  reexpansion  for  the  transient  parameters  exhibit  the 


essential  features 
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Figure  16.  Transient  currents  for  a plane  wave  with  a Gaussian  tine 
history  at  0 degrees  Incidence.  Solid  lines  are  from  TWTD,  dashed  lines 
are  fron  the  expansion  of  SEM  parameters  obtained  from  an  integral  equa- 
tion technique,  and  the  dot-dashed  lines  are  an  expansion  of  the  SEM 
parameters  derived  from  transient  data. 


Observation  zone  12 


Figure  17.  Transient  currents  for  a plane  wave  with  a Gaussian  time 
history  at  30  degrees  incidence.  Solid  lines  are  from  TWTD,  dashed 
lines  are  from  the  expansion  of  SEM  parameters  obtained  from  an  Integral 
equation  technique,  and  the  dot-dashed  lines  are  an  expansion  of  the 
SEM  parameters  derived  from  transient  data. 
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Figure  18.  Transient  currents  for  a plane  wave  with  a Gaussian  tine 
history  at  60  degrees  incidence.  Solid  lines  are  from  TWTD,  dashed 
lines  are  from  the  expansion  of  SEM  parameters  obtained  from  an  integral 
equation  technique,  and  the  dot-dashed  lines  are  an  expansion  of  the  SEM 
parameters  derived  from  transient  data. 
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4.3  Reexpansion  for  Double  Exponential  Excitation 

Using  a plane  wave  with  a double  exponential  time  history  as  the 
excitation,  the  derived  SEM  data  base  was  reexpanded  for  0,  30,  and  60 
degrees  incidence.  The  waveform  pole  contributions  are  included  in  these 
expansions.  Thus  we  may  compare  early  time  results  as  well  as  the  late 
time.  The  current  was  observed  at  observation  zones  12,  18,  and  25  for 
each  excitation.  These  results  were  compared  with  a reexpansion  in  terms 
of  SEM  parameters  obtained  by  a method  of  moments  computation. 

The  transient  current  observed  at  the  three  observation  segments  at 
0 degrees  incidence  is  shown  in  Figure  19.  We  can  see  from  this  Figure 
the  transient  extracted  reexpansion  and  the  reexpansion  for  the  method  of 
moments  data  compare  very  well,  in  both  the  early  and  the  late  time. 

Figure  20  displays  the  transient  current  at  30  degrees  for  the  three 
observations.  Here  we  see  good  agreement  in  the  late  time  for  observation 
zones  18  and  25,  but  some  minor  departure  at  zone  12  for  both  the  early 
and  late  times.  These  departures  can  be  attributed  to  the  imperfect 
extracted  modes  and  normalization  constants. 

The  three  observation  zones  for  60  degrees  are  shown  in  Figure  21. 
Here  we  see  good  agreement  among  the  two  cases  for  late  time  with  some 
departure  among  the  two  in  the  early  time. 

In  summary  the  reexpansion  of  the  transient  extracted  data,  in  gen- 
eral, compares  very  well  with  the  transient  currents  computed  by  TWTD, 
and  the  expanded  SEM  data  derived  by  the  method  of  moments.  We  therefore 
conclude  that  extrapolation  of  transient  response  data  to  new  excitations 
using  the  Prony-based  SEM  extraction  procedure  is  feasible. 
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Figure  19.  Transient  currents  for  n plane  wove  wltli  n double  exponcnHal 
tine  history  at  0 degrees  incidence.  Solid  lines  are  from  an  expansion 
of  the  SEM  parameters  derived  from  transient  data,  and  the  dashed  lines 
are  from  an  expansion  of  the  SEM  parameters  obtained  from  an  integral 
equation  technique. 


Figure  20.  Transient  currents  for  a plane  wave  with  a double  exponential 
time  history  at  30  degrees  incidence.  Solid  lines  are  from  an  expansion 
of  the  SEM  parameters  derived  from  transient  data,  and  the  dashed  1 incs 
are  from  an  expansion  of  the  SOI  parameters  obtained  from  an  integral 
equation  technique. 


Figure  21.  Transient  currents  for  a plane  wave  with  a double  exponential 
time  history  at  60  degrees  incidence.  Solid  lines  are  from  an  expansion 
of  the  SEM  parameters  derived  from  transient  data,  and  the  dashed  lines 
are  from  an  expansion  of  the  SEM  parameters  obtained  from  an  integral 
equation  technique. 


CHAPTER  V 


CONCLUSIONS 

Using  the  method  presented  in  this  document  it  appears  feasible  to 
extract  the  SEM  description  from  spatially-sampled  transient  data.  To 
obtain  a complete  SEM  description  for  a system  several  data  completeness 
considerations  must  be  taken  into  account.  They  are  as  follows: 

- For  extraction  of  the  SEM  description,  the  magnitude  of  the  spectrum 
of  the  forcing  function  at  the  highest  frequency  of  interest  must  be 
well  above  the  noise  level  to  insure  this  frequency  gets  excited; 

- The  spatial  placement  of  the  generator  must  be  such  that  all  modes 
of  interest  are  excited;  and 

- The  extrapolation  of  the  SEM  description  cannot  be  carried  out  if 
the  expansion  extends  beyond  the  band  limitations  of  the  descrip- 
tion, or  if  the  excitation  couples  data  which  were  not  coupled  in 
the  data  from  which  the  description  is  derived. 

The  quality  of  data  extracted  is  directly  dependent  on  the  quality 
of  the  transient  responses.  For  center -excitation  conditions  the  coupling 
is  a maximum  for  even  modes,  and  the  resulting  SEM  description  was  very 
favorable.  For  off-center  excitation  the  coupling  is  not  as  strong,  there- 
fore  the  results  were  weaker  than  for  a center  excitation.  In  general,  we 
observe  that  the  information  content  of  the  original  data  is  transformable 
to  the  SEM  form. 

Results  of  the  extrapolation  to  different  excitations  were  favorable. 
This  proved  true  even  when  clearly  discernable  errors  were  present  in  the 
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